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As discussed in the accompanying article yj], Casimir first predicted that zero-point quantum 
fluctuations of an electromagnetic (EM) field give rise to an attractive force between two closely 
spaced perfect conductors M\. In this idealized situation, a macroscopic physical manifestation, 
i.e. a force, arises purely from the EM quantum fluctuations. Similarly, one would expect an 
analogous force to arise from quantum fluctuations in a quasiparticle vacuum in a zero-temperature 
dilute Bose-Einstein condensate (BEC). 

In this article, an argument is presented for the existence of a Casimir-like drag force on a 
stationary object in a slow-moving BEC due to the scattering of quantum fluctuations. This force 
would be the dominant effect at speeds well below that determined by the Landau criterion Q3[| 
(the derivation of which ignores the possibility of scattered quantum fluctuations). No direct EM 
analogy can be drawn for this situation because of the Lorentz invariance of the EM vacuum, i.e. 
no absolute rest frame (where the relative motion of an object can be measured) exists for an EM 
vacuum. Before exploring the Casimir-like drag force, the connection between the familiar EM 
Casimir force and the Casimir-like force in a static BEC situation shall first be made explicit. 

The following estimate of a Casimir-like pressure pc in a weakly interacting BEC will serve to 
motivate the existence of a force in a BEC vacuum due to quantum fluctuations. Consider two in- 
finitely thin and infinitely repulsive plates (acting in an analogous role to perfect conductors in an 
EM field) spaced a distance b apart in the ^-direction, 111 ill immersed in a three-dimensional con- 
densate characterized by an interparticle contact pseudopotential, gfi( 3 \r), where g is determined 
by the 2-particle positive scattering length a and the mass m of the atoms, i.e. g = 4:Tch 2 a/m. 



Let us assume the condensate is dilute such that a/ po« 3 "C 1 where p is the condensate density. 
In the homogeneous dilute condensate (without confining walls), the ground state energy can be 
written Ji li as E g (k) = p 2 gV/2 + l/2£ (k#0) [£ B (fc) - E (k) - Po g] [12] where k = |k|, 
E (k) = h 2 k 2 /(2m), E B (k) = [(E (k)) 2 + 2p gE (k)] 1/2 is the Bogoliubov dispersion relation, 
po is the condensate density, and V is the volume of the condensate. The pressure generated by the 
leading order interaction term is given by p = gpl/2. Subsequent terms describe the contribution 
from quantum fluctuations, or quantum depletion, to the ground state energy, which gives rise to a 
Casimir-like force. Assuming the distance between the confining walls to be much greater than the 
healing length, i.e. b ^> k^ 1 = (87rp a) _1 ^ 2 1130, Casimir's original calculation can be followed 
to determine the change in energy per area between the homogeneous case, which displays a con- 
tinuous momentum spectrum, and the case with confining walls, which has a discrete momentum 
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spectrum along z. The derivative with respect to the distance yields the Casimir-like pressure 

d ( 1 °° f°° \ 

Pc = --^ laxity + X{n) - J dnx{n)\ , (1) 

where x{ n ) = V 47r f k r dk r [E B (n) - E (n) - p g], k 2 (n) = k 2 . + {nn/b) 2 , and A kr is the 
ultraviolet cutoff where the simple assumption of a delta-function contact pseudopotential breaks 
down. Using the Euler-MacLaurin formula and taking A k — > 00 gives the Casimir-like pressure in 
a dilute condensate to leading order as 

where the derivative is with respect to n and the speed of sound in this system is given by c s = 



y/po9/m. 

The connection with the Casimir force is clear: this expression for the Casimir-like pressure is 
exactly one half of the original Casimir pressure arising from the zero-point fluctuations of the EM 
vacuum, with the speed of sound replacing the speed of light. This connection can be understood 
in light of the observation that the lower wavenumber modes dominate the Casimir-like pressure 
(the higher modes cancel out as the discrete momentum spectrum becomes more like a continuous 
spectrum). At long wavelengths the ground state energy is approximately E g (k) = p^gV/2 + 
l/2^(k^o)[^ c s^ — Po5'] 5 which, aside from an unimportant constant factor and the substitution 
of the speed of light with c s , is the same harmonic oscillator energy spectrum that describes the 
EM vacuum. The factor of one half between the BEC case and the EM case arises as there are 
two transverse polarization states for each EM mode whereas there is only one polarization state 
for the phonon modes in a BEC. Although for a dilute condensate this force is 10 6 times smaller 
than the Casimir force in an EM vacuum and finite-size effects present in experiments have not 
been included, this force could in principle be measured, providing direct evidence of the quantum 
fluctuations in these systems. 

The next step is to show that a drag force also arises from these fluctuations in a slow-moving 
BEC flow around a stationary object. Since this force arises from quantum fluctuations, one needs 
to go beyond the usual mean field picture and treat the fluctuations as scattered waves. An outline 
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of this calculation for a general object is presented below (for more details see [8, 

In general, a force on a stationary object described by a potential V(r) (r represents a vector 
quantity throughout this paper) can be written in second quantized notation at zero temperature as 

p = - I d 3 r(^(r)[VV(r)]4j(r)) T=0 , (3) 



where %jj(r) and ^(r) are field operators that describe the weakly interacting BEC flow and obey 
the standard boson commutation relations (the expectation value at T = is defined below). 
Landau, based on his phenomenological theory that did not explicitly consider the quantum fluc- 
tuations, showed the moving ground state to decay by the emission of quasiparticles only above 
a critical speed; below that critical speed the flow would be dissipationless [3]. Consequently, 
according to Landau's theory a stationary object in a superfluid flowing below this critical velocity 
would remain metastable and not experience any force. 

Landau's argument, however, only directly applies to a uniform condensate (containing in- 
finitesimally small test particles that do not perturb the condensate flow) whose state is derived 
from the immobile ground state by a Galilean transformation. The critical velocity in a non- 
uniform system is more complicated as the quantum state is not derived by a simple uniform 
Galilean boost. At zero temperature it can generally be assumed that the condensate, whether uni- 
form or non-uniform, is well described by a macroscopic classical field, i.e. ip(r) = ^^(r). The 
force on a stationary object in a moving flow is then simply given by 

F GPE = - J d 3 r\¥°\r)\ 2 VV(r) (4) 

and the behavior of this field in a time-independent situation is given by the Gross-Pitaevskii 
equation (GPE) 

(f + V(r) - fi)¥°\r) + g\¥ 0) (r)\ 2 ¥°\r) = 0, (5) 

where [i is the chemical potential, T = — + ifrc-^ + ^mc 2 and c is the speed of the flow far 
from the object. In the non-uniform case, it is the maximum local velocity that should be com- 
pared to Landau's critical velocity. Using the GPE, one can show that Fgpe is zero only for a 
flow velocity at infinity that is associated with a maximum local flow velocity that does not exceed 
Landau's critical velocity. Although the critical velocity, as measured at infinity, is modified in the 
non-uniform case, the principle of Landau's argument still applies. However, if quantum fluctua- 
tions — beyond the scope of the GPE — are included then even such a modified critical velocity 
defining a transition to frictionless flow no longer necessarily applies to this situation. Hereafter, to 
avoid confusion in the discussion below, no explicit distinction will be made between the modified 
critical velocity in the non-uniform case and Landau's critical velocity in the homogeneous case. 

Below, the calculation for the force due to the quantum fluctuations at T = is described. The 
assumption is made that the superfluid is well approximated by a weakly interacting Bose-Einstein 
condensate. The Bose field operator that describes the flow is approximated as a large macroscopic 
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classical field (as above) plus a small quantum fluctuation operator (the small parameter in this case 
is given by the diluteness parameter ^/ poa 3 [6]), i.e. 

^S(r) = ^«(r) + (j>(r), (6) 

where is the macroscopic condensate that has been modified by the quantum fluctuations. 
The force acting on a stationary potential can then be written as 

d 3 r(\¥ l \r)\ 2 + (ft(r)j>(r)) T=0 )VV(r). (7) 

Next, cf> (5|] is expanded in terms of quasiparticle operators a k and a\ that, by definition, obey 
the standard boson commutation relations [a k , a\,\ = 8 kjk > where 8 kjk i is the Kronecker delta 
function. The quasiparticle operators are weighted by the quasiparticle amplitudes u k {r) and v k (r) , 
i.e. 

<K r ) = ( u k( r )®k - vl(r)a£) , (8) 

(k?40) 

where the sum is taken over all excited states and excludes the condensate mode and * denotes the 
complex conjugation. The assumption is made that the quasiparticle amplitudes satisfy the two 
sets of differential equations [6], also known as the Bogoliubov equations, 

tu k {r) - g(*®) 2 v k (r) = E k u k (r) (9) 

t*v k (r) - g(¥°>) 2 u k (r) = -E k v k (r), (10) 

where * denotes the complex conjugation, E k is the eigenvalue associated with momentum state 
k, and C = T + V(r) — \i + 2g\^/^\ 2 . The quasiparticle operators would then diagonalize the 
weakly interacting Hamiltonian. In other words, weakly interacting particles are transformed into 
non-interacting quasiparticles. 

The behavior of the condensate modified by the quantum fluctuations at T = is given by the 
generalized GPE QVI] 



(f+V(r)-fi)¥ 1 \r)+g\¥ 1 \r)\ 2 ¥ 1 \r)+ ^ 2g\v k (r)\ 2 ¥ l \r)-g ^ u k (r)vl{r)m {1 \r)-f(r)m {1 \r) 
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The term proportional to ^( k ^ ) u k( r ) v k ( r ) 1S ultraviolet divergent because of the contact poten- 
tial approximation and must be renormalized. The last term, f(r)^/^(r), ensures the orthogonal- 
ity between the excited modes and the condensate and is given by /(r) = 5^(k^o) c k v k ( r )> where 
c k = gj d 3 r|#«(r)| 2 (#«*(rK(r) + V®(r)v k (r)). 



Since the fluctuations are assumed to be small, the fluctuation part of the condensate can be sep- 
arated from the macroscopic part described by the GPE, i.e. ^W(r) = \E^ )(r) + J^(k^o) ^k(r), 
where 8^f.{ r ) is associated with a given mode k. Therefore, since it is known that F GPE = at 
speeds lower than the critical velocity, at these low speeds the force (at the dominant order) on a 
stationary object can be written as 

Ffi uc = ~J2[ d 3 r^ i0 \r)5^ k (rY + c.c. + \v k (r)\ 2 )VV(r), (12) 
(k^o) J 

which is proportional to the diluteness parameter a/ poa 3 . 

To complete the calculation of the Casimir-like force, first the coupled Bogoliubov equations 
(eqns [9] [101) must be solved for the quasiparticle amplitudes, then the generalized GPE must be 
solved for 5^k(r), and finally all excited states must be integrated over. The Bogoliubov equa- 
tions can be treated as a scattering problem in which the only incoming waves considered are the 
incident homogeneous fluctuations (incoming scattered waves are excluded); exponentially grow- 
ing scattered waves are excluded. Temperature is not well defined in this scattering problem, so 
instead one can impose in terms of measurable quantities (specifically the spatio-temporal correla- 
tion functions) a condition far from the potential that determines the quantum state in this problem 
(see [8|] for details) such that the state is annihilated by a. Therefore, it remains convenient to 
denote that state by T = 0'. The expectation values written in terms of quasiparticle amplitudes 
thus take the form (ft(r)4>(r)) T=0 = Xl(k^o) \ v k{ r )\ 2 , for example. 

The calculation in the Born approximation is described in [8] for a specific potential that varies 
only in the flow direction, and in |90 for a 3-D weak point impurity, both in an infinitely extended 
medium. In the latter case, the force is shown to be linearly proportional to the flow velocity. A 
force is found to exist at all flow velocities in both cases. 

Ffiuc, a direct physical manifestation of the zero-point fluctuations, dominates at velocities 
below Landau's critical velocity because, due to the Landau criterion, there is no contribution 
described by the GPE (eq. [5]) to the overall force. Ff iuc comes not from a change in the number 
of quasiparticles, but from a change in the associated amplitudes resulting from the scattering 



condition necessary to solve the Bogoliubov equations. Assuming this force to be dissipative [14], 
energy is transferred by the scattering of fluctuations (analogous to wave drag in classical fluid 
mechanics) rather than by the creation of quasiparticles at higher flow velocities or the scattering 
of quasiparticles at non-zero temperature. 

Although the existence of a drag force on a localized object in an infinitely extended BEC 



flow seems to imply that persistent currents cannot exist, this effect can be consistent with the 
observed persistent currents, where superfluids placed in a closed loop are able to flow practically 
endlessly without any detectable friction. Consider a closed loop of flowing dilute condensate 
at zero temperature contained by perfectly smooth walls (we consider rough walls below) with 
an impurity fixed in the flow. From the discussion above, one would expect (ignoring typical 
transient effects, unrelated to quantum fluctuations, that arise from the act of placing the impurity 
in the flow) that for a short time interval from the moment the impurity is introduced in the flow, 
this system would act as if the medium were infinitely extended, i.e. the impurity would feel a 
force and the flow would dissipate energy. However, the scattering fluctuations experience the 
finitude of the system at time scales on the order of the characteristic length of the system divided 
by the speed of sound. The scattering fluctuations then re-interact with the impurity, eventually 
cancelling out any drag force. A certain amount of scattered fluctuations must be present before 
this force can cancel out. If one includes higher order terms in the Hamiltonian that describe the 
quasiparticle interactions (i.e. the Beliaev and Landau damping terms), one can assume the system 
will relax to a local thermal equilibrium. In other words, the scattered fluctuations will relax to a 
small amount of normal gas. Therefore, one would expect there to be two relevant observables — 
a new time scale for the finitude of the system to be felt and a small amount of normal gas that 
must be produced to counteract the force on the fixed impurity far from the object. 

Finally, one expects the scattering of quantum fluctuations off a rough surface to have dramatic 
effects on the superfluid boundary conditions. As discussed in [Ji]], in addition to the total normal 
mass flux condition for a solid wall, n- (j n +js) = where n is the normal vector to the surface and 
j n and j s are the mass fluxes of the normal and superfluid, the scattering of quantum fluctuations 
off a rough surface leads to the following non-equilibrium hydrodyamic boundary condition 

n-(J n - 1) = a(T b - T) + (3v 2 sl (13) 

where T b is the temperature at the boundary, and T is the temperature of the fluid nearby. This 
boundary condition holds to leading order and has important predictions, such as a new critical 
speed and that there should be a temperature difference between a moving superfluid and the 
boundary, which implies that there should always be at least a small amount of normal fluid de- 
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tectable in immersed torsion pendulum experiments. 
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